"Reservoir model" for shallow modulation-doped digital magnetic quantum wells 



Henrique J. P. Freire^'l^] and J. Carlos Egues^'|] 

^ Departamento de Fisica e Informdtica, Instituto de Fisica de Sao Carlos, 
Universidade de Sao Paulo, Caixa Postal 369, 13560-970 Sao Carlos-SP, Brazil 
'^Department of Physics and Astronomy, University of Basel, 
Klingelbergstrasse 82, CH-4056 Basel, Switzerland 
(Dated: July 13, 2002) 

Digital Magnetic Heterostructures (DMH) are semiconductor structures with magnetic mono- 
layers. Here we study electronic and magneto-transport properties of shallow modulation-doped 
(ZnSe/ZnCdSe) DMHs with spin-5/2 Mn impurities. We compare the "reservoir" model, possibly 
relevant to shallow geometries, to the usual "constant-density" model. Our results are obtained 
by solving the Kohn-Sham equations within the Local Spin Density Approximation (LSDA). In 
the presence of a magnetic field, we show that both models exhibit characteristic behaviors for 
the electronic structure, two-dimensional carrier density, Fermi level and transport properties. 
Our results illustrate the relevance of exchange and correlation effects in the study shallow 
heterostructures of the group II- VI. 
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Digital Magnetic Heterostructures (DMHs) are state 
of the art layered semiconductor structures in which 
(quasi-) two-dimensional distributions of magnetic mo- 
ments are restricted to equidistant planes within a quan- 
tum well ||l|, ||. Spin-dependent quantum effects are 
pronounced in these systems. For example, in an ex- 
ternal magnetic field the s-d exchange interaction ||] 
between itinerant electrons and those of the magnetic 
impurities is responsible for a giant Zeeman effect up 
to two orders of magnitude larger than the ordinary 
one. These large energy splittings are readily ob- 
served by magneto-photoluminescence [0, In ad- 
dition, the successful achievement of high-doping car- 
rier densities in Q has enabled magnetotransport and 
magneto-photoluminescence measurements in these sys- 
terns i I, |. 

In this work we investigate spin-dependent properties 
of shallow DMHs, Fig. 1(a). "Shallow" here implies that 
the confining potential is weak enough, e.g. ~ 25 meV 
for a ~ 10.5 nm wide well [Fig. 1(b)], so as to have only 
a few confined subbands. For such geometries, the adja- 
cent n-doped regions with densities ~ 10^^ cm~^ provide 
enough carriers to fill up all confined levels in the well, 
thus serving as electron reservoirs. We describe these 
shallow DMHs using a "reservoir model" with a constant 
Fermi level pinned to the chemical potential of the n- 
doped regions |l^. Such a model has also been ap- 
]^d to GaAs/GaAlAs [|l], [l2[|l|l and to ZnSe/ZnCdSe 
Q systems. We calculate the magnetic-field depen- 
dent subband structure of our shallow DMHs by using 
Density Functional Theory within a Local Spin Den- 
sity Approximation (DFT/LSDA). As recently shown. 



exchange-correlation effects are important in shallow II- 
VI wells We determine Landau-level fan diagrams, 
two-dimensional electron densities, and in-plane trans- 
verse resistivities. 
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FIG. 1: (a) Layered structure of a digital magnetic quan- 
tum well heterostructure. Lateral n-doped regions (impurity 
concentration of 1.2 x lO^'^ cm~'^) provide electrons that fill 
the confined energy levels. Self-consistent potential profile for 
majority spin down electrons (at B = 1 T) are shown for the 
"reservoir" (b) and "constant density" (c) regimes. Each case 
has only one confined subband bellow the Fermi level. In (b) 
total electronic density n(z) penetrates into the lateral bar- 
riers allowing exchange and correlation to contribute to the 
potential profile far from the well center. 
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FIG. 2: Magnetic field dependence of calculated Landau levels, two-dimensional densities, and transverse resistivities for both 
the reservoir (upper panels) and the constant density (lower panels) models. Within the reservoir model the Landau levels 
(a) are pinned to the Fermi level ef of the reservoirs, thus originating sawtooth-profile LL fan diagrams, ef (b) is pinned to 
the highest occupied LL in the constant-density model and oscillates as the magnetic field is increased. The s-d enhanced LL 
spin splitting gives rise to the spin polarization of the 2DEG at low B fields [(c) and (d)]. Note the oscillatory behavior of the 
two-dimensional density within the reservoir model (c). The transverse resisitivities p^x display Shubnikov-de Haas oscillations 
in both regimes (e) and (f); however, plateaus pxx = are seen only in (e) since we do not account for localized states in this 
model. 



We compare our results within the "reservoir model" 
with those obtained with the conventional "constant- 
density" model, in which deep enough wells can confine 
all available carriers. In this case the n-doped regions 
adjacent to the well are fully depleted. We configure the 
system geometry for each case by suitably choosing the 
height of the confining barriers at the ZnSe/ (Zn,Cd)Se in- 
terfaces: 70 meV and 25 meV for the "constant density" 
and the "reservoir" models, respectively. The idea is to 
identify relevant contrasting features in these models to 
eventually better describe real shallow DMHs. 

We calculate the magnetic-field dependent subband 
structure of our system within Density Functional The- 
ory (DFT). We solve the Kohn-Sham (KS) equations in 
the Local Spin Density Approximation (LSDA) using the 
parametrization of VWN . We compare our LSDA re- 
sults to those obtained within the Hartree approximation 
to assess the role of exchange and correlation (XC) effects 
in the properties of these novel systems. We consider the 
system at zero temperature and perform the calculations 
in the effective-mass approximation framework. We also 
make the usual assumption that we can interpret the re- 
sulting KS eigenvalues as subband energies of the system. 
Consequently, the total energy of each Landau level (LL) 
is the sum of the subband energy Si^a-^ and the Landau 
cyclotron energy {n + \/2)hwc 



£i,n,o, [B) = £i,a, {B) + [n + - \nujc + -gusazB. (1) 



In (|l|), the subband energy ei^cr^{B) corresponds to the 
longitudinal confined electronic motion and the second 
term to the quantized transversal motion due to the mag- 
netic field. For completeness we also add the Zeeman en- 
ergy, although its contribution is very small when com- 
pared to the s-d exchange interaction. 

Panels (b) and (c) of Fig. 1 show the resulting self- 
consistent potential profiles for both "reservoir" (a) and 
"constant-density" (b) regimes. As the magnetic field in- 
creases, the self-consistent potential and the correspond- 
ing electronic structure change due to i) magnetic-field 
dependent contributions to the confining potential and 
a) the rearrangement of electrons in the confined lev- 
els (magnetic-field dependent Landau-level degeneracy 
eB/h). 

Figures 2(a) and (b) show the Landau- level fan dia- 
grams for both regimes and the magnetic field depen- 
dence of the Fermi energy £^7-. In the "reservoir" regime 
case the Fermi level is constant at the chemical potential 
of the reservoir (n-doped regions) while in the "constant 
density" case ep oscillates as a function of B. In the 
first case, all confined electronic levels are fully occupied 
because n-doped regions can always provide enough elec- 
trons to the QW, thus being only partially depleted. For 
this reason the LLs are pinned to the Fermi level (see 
plateaus at £f)- In the "constant density" case the n- 
doped regions are completely depicted but the amount 
of electrons provided to the QW is not sufficient to fill 
all electronic levels. Consequently, the Fermi energy is 
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pinned to the highest occupied LL, and oscillates as a 
function of the magnetic field [Fig. 2(b)]. Note that in 
this case is pinned to the LL, not the contrary as in the 
"reservoir" regime. This is why the curves in former case 
exhibit a sawtooth profile not present in the latter. Note 
also that in the "reservoir" ( "constant density" ) case the 
two-dimensional density n2D changes (is constant) as a 
function of B [Fig. 2(c) and (d)]. 

When compared to pure Hartree calculations (c.f. Ref. 
[^) we note that XC plays an important role in such 
group II- VI heterostructures. Within LSDA calculations 
the LL spin splitting is enhanced, the two-dimensional 
density is increased, and the magnetic field at which the 
2DEG becomes fully polarized is decreased. 

We also calculate the magnetic field dependence of the 
in-plane resistivities (perpendicular to the magnetic field 
and growth axis). We use Ando's model for the lon- 
gitudinal conductivity axx ^"^^ Drude's model for 
the transversal component of the conductivity, Uxy — 
en2D/B. These simple models enable us to illustrate 
general magnetotransport features within these n-doped 
magnetic structures and to compare the two regimes for 
the density. In Fig. 2(e) and (f) we plot our results for 
the in-plane longitudinal resistivity (p = f ^^) which is a 
relevant experimental quantity. 

Both "reservoir" and "constant density" regimes dis- 
play the usual oscillations in p^x (Shubnikov-de Haas os- 
cillations), but the former has regions of zero resistivity 
[Fig. 2(e)]. In the "reservoir" regime the QW is always 



filled with the maximum number of electrons possible, 
which is defined by both the magnetic-field dependent 
density of states and the reservoir chemical potential. 
Within this self-consistent equilibrium, there are mag- 
netic field ranges at which the center of the LLs are not 
on the Fermi level [Fig. 2(a)], and consequently, the re- 
sistivity is minimum [Fig. 2(e)]. In the "constant den- 
sity" case, on the other hand, the Fermi level is always 
pinned to highest occupied LL and the minimum resis- 
tivity dips [Fig. 2(f)] occur at magnetic fields that ep 
jumps between LLs [Fig. 2(b)J. At those magnetic fields, 
ep lies on the tail of a LL l|l8| so that there is only a 
small number of conducting electrons and, consequently, 
a correspondingly small resistivity. In passing, we note 
that in the standard model for the integer quantum hall 
effect (IQHE) localized states in broadened LLs 

are responsible for the Fermi-level pinning |2^; the phe- 
nomenological inclusion of these states (gaussian density 
of states) gives rise to plateaus pxx = 0. 

Finally, we should mention that we are currently in- 
vestigating the temperature dependence of the transverse 
magnetoresistivity |2^; a more detailed analysis, beyond 
the scope of this communication, will be addressed in a 
future publication. This additional study should allow 
us to contrast the two models discussed here in terms of 
agreement with experimental data. 
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